Aft~oiQ-feG6 


Synchronous  Averaging  of  Helicopter 
Tail  Rotor  Gearbox  Vibration:  Phase 
Reference  Considerations 


D.M.  Blunt 
DSTO-TR-0739 


APPROVED  FOR  PUBLIC  RELEASE 

©  Commonwealth  of  Australia 

Department 


♦ 


O  F 


Defence 


DEFENCE  SCIENCE  AND  TECHNOLOGY  ORGANISATION 


Synchronous  Averaging  of  Helicopter  Tail  Rotor 
Gearbox  Vibration:  Phase  Reference  Considerations 


DM.  Blunt 

Airframes  and  Engines  Division 
Aeronautical  and  Maritime  Research  Laboratory 

DSTO-TR-0739 


ABSTRACT 

Synchronous  averaging  requires  an  accurate  phase  reference  (tachometer)  signal.  In 
helicopter  transmissions,  such  a  signal  can  usually  be  obtained  or  derived  horn  an 
engine  or  main  rotor  gearbox  accessory,  but  is  not  as  readily  available  for  the  tail  rotor 
gearbox.  This  report  examines  whether  a  separate  tail  rotor  gearbox  phase  reference 
signal  is  necessary  by  investigating  the  relative  jitter,  due  to  dynamic  tail  drive  shaft 
twist,  between  phase  reference  signals  obtained  from  the  main  and  tail  rotor  gearboxes 
in  a  S-70A-9  Black  Hawk  helicopter. 
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Synchronous  Averaging  of  Helicopter  Tail 
Rotor  Gearbox  Vibration:  Phase  Reference 

Considerations 


Executive  Summary 

Synchronous  averaging  is  a  signal  processing  technique  that  significantly  enhances  the 
ability  to  detect  and  diagnose  faults  in  rotating  machinery  by  'extracting'  the  vibration 
for  a  particular  shaft  from  the  overall  gearbox  vibration.  It  is  based  on  the  principle 
that  the  forced  vibration  of  a  shaft,  and  the  components  attached  to  and  rotating  with 
it,  will  be  periodic  with  the  shaft  rotation.  The  technique  requires  an  accurate  phase 
reference  to  synchronise  the  sampling  of  the  vibration  signal  with  the  rotational 
frequency  of  the  shaft.  The  synchronous  samples  are  then  divided  into  segments  of  one 
rotational  period,  and  ensemble-averaged  to  attenuate  the  non-synchronous 
frequencies. 

Because  it  is  often  difficult  to  obtain  a  phase  reference  signal  directly  from  a  helicopter 
tail  rotor  gearbox,  this  investigation  was  conducted  to  determine  whether  such  a  direct 
signal  was  necessary,  or  if  an  indirect  phase  reference  signal  from  the  main  rotor 
gearbox  would  be  sufficient.  The  investigation  focused  on  the  possibilty  for  dynamic 
twisting  of  the  long  tail  drive  shaft,  and  how  this  might  produce  phase  errors  at  the  tail 
rotor  shaft.  These  errors,  even  if  very  small,  have  the  potential  to  significantly  reduce 
the  sensitivity  of  fault  detection  methods  based  on  the  synchronous  averaging  process. 

The  data  from  a  recent  flight  trial  in  an  Australian  Army  S-70A-9  Black  Hawk 
helicopter  were  analysed  to  determine: 

a)  the  relative  differences  between  main  and  tail  phase  reference  signals,  and 

b)  the  effect  of  these  differences  on  the  synchronous  averages  of  the  tail  gearbox 
output  shaft  vibration. 

It  was  found  that  during  forward  flight  the  dynamic  twisting  of  the  tail  drive  shaft 
appears  to  be  minimised,  but  that  during  hover,  and  particularly  tail  rotor  turns,  there 
appears  to  be  a  significant  amount  of  twisting.  This  was  found  to  produce  a 
consequential  degradation  of  the  synchronous  averages. 

It  was  also  discovered  that,  in  the  S-70A-9,  forward  flight  induces  undesirable 
vibration  impulses  into  the  tail  gearbox  vibration  signal.  These  have  subsequently  been 
confirmed  to  exist  in  the  S-70B-2  as  well.  The  cause  of  the  impulses  is  speculated  to  be 
the  aerodynamic  interaction  of  the  main  rotor  tip  vortices  on  the  tail  rotor  blades.  The 
impulses  are  undesirable  from  a  vibration  analysis  point  of  view  because  of  their 
potential  to  mask  faults  within  the  gearbox. 


It  is  therefore  recommended  that: 

a)  a  direct  phase  reference  signal  from  the  tail  rotor  gearbox  should  always  be  used, 
and 

b)  S-70/H-60  series  tail  rotor  gearbox  vibration  should  only  be  analysed  during 
hover. 

These  recommendations  have  a  direct  outcome  for  the  AMRL-developed  transmission 
vibration  monitoring  system  currently  being  implemented  in  all  Royal  Australian 
Navy  S-70B-2  Seahawk  and  SK-50  Sea  King  helicopters.  This  system  consists  of 
vibration  transducers  hard-wired  into  the  aircraft,  and  a  carry-on/ carry-off  data 
acquisition/ analysis  unit.  Tail  rotor  photocells  are  being  fitted  to  both  types  of 
helicopters,  and  tail  rotor  gearbox  vibration  data  will  only  be  acquired  during  hover 
out  of  ground  effect. 
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1.  Introduction 

Synchronous  averaging  is  a  technique  regularly  employed  in  the  analysis  of  helicopter 
transmission  vibration  to  provide  estimates  of  the  vibration  from  individual  shafts  and 
their  associated  components.  It  is  fundamentally  based  on  the  principle  that  the  forced 
vibration  of  a  shaft,  and  the  components  fixed  to  and  rotating  with  it  (eg  gears),  will  be 
periodic  with  the  shaft  rotation.1  Synchronous  averaging  requires  an  accurate  shaft 
phase  reference  (tachometer)  signal  to  divide  the  gearbox  vibration  signal  into  periods 
of  exactly  one  shaft  revolution.  These  segments  are  then  ensemble-averaged  to 
attenuate  the  non-synchronous  components.  The  resultant  average  is  treated  as  if  it 
were  the  time-domain  vibration  signal  for  one  revolution  of  the  isolated  shaft  with  its 
attached  components. 

As  would  be  expected,  synchronous  averaging  is  very  dependent  on  the  accuracy  of 
the  phase  reference  signal.  Any  phase  jitter,  or  error,  in  this  signal  will  make  the 
synchronous  vibration  appear  slightly  non-synchronous,  and  cause  it  to  be  attenuated 
during  the  averaging  process. 

Since  gearbox  shafts  are  kinematically  linked  by  their  gear  ratios,  a  single  phase 
reference  signal  will  usually  be  sufficient  to  synchronously  average  the  vibration  of  all 
the  gearbox  shafts.  The  use  of  this  signal  may  also  be  extended  to  external  drivetrain 
components,  but  the  further  the  components  are  from  the  location  of  the  phase 
reference,  the  more  likely  it  is  for  the  flexibility  of  the  drive  system  to  introduce 
additional  phase  errors.  This  will  be  even  more  likely  if  the  drivetrain  is  subject  to 
torque  fluctuations.  These  additional  phase  errors  will  be  superimposed  on  top  of  those 
inherent  in  the  phase  reference  signal,  and  degrade  the  synchronous  average  even 
further. 

In  helicopter  transmissions,  the  gearboxes  that  will  be  most  susceptible  to  this 
phenomenon  are  those  at  the  end  of  the  long  tail  drive  shaft:  that  is,  the  tail  rotor 
gearbox  and,  if  it  exists,  the  intermediate  gearbox  at  the  base  of  the  tail  pylon.  This  is 
because  phase  reference  signals  can  usually  be  obtained,  or  derived,  from  an  engine  or 
main  rotor  gearbox  accessory,  such  as  an  AC  generator,  but  are  not  readily  available 
for  the  intermediate  or  tail  rotor  gearboxes.  It  should  also  be  noted  that,  in  most 
helicopters,  the  tail  drive  shaft  is  typically  made  of  three  or  four  shorter  shafts.  These 
are  supported  by  hanger  bearings  mounted  in  soft  viscous  dampers,  and  connected  by 
flexible  couplings,  both  of  which  may  contribute  to  the  overall  flexibility  of  the  shaft. 

To  investigate  this  potential  problem,  a  photocell  was  fitted  to  the  tail  rotor  of  an 
Australian  Army  S-70A-9  Black  Hawk  helicopter  during  a  recent  trial  installation  of  a 
transmission  vibration  monitoring  system  [1].  This  supplemented  a  phase  reference 
signal  derived  from  an  AC  generator  geared  into  the  main  rotor  transmission.  Data 


1  This  excludes  rolling  element  bearings  because  of  slip. 
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were  captured  under  a  number  of  different  flight  conditions,  and  the  relative  phase 
jitter  between  the  two  signals  was  examined  by: 

a)  calculating  the  difference  in  tail  rotor  phase  angle  derived  from  each  signal,  and 

b)  comparing  the  synchronous  averages  of  the  tail  rotor  gearbox  output  shaft 
computed  using  each  signal. 


2.  Theory 

A  brief  description  of  synchronous  averaging  is  presented  below.  For  a  more  detailed 
examination,  particularly  on  the  effects  of  various  interpolation  methods  on  the  digital 
re-sampling  of  vibration  data,  the  reader  is  referred  to  Forrester  [2].  This  description  is 
followed  by  a  simple  mathematical  model  of  the  attenuation  effect  of  reference  phase 
jitter  on  a  synchronous  average. 

2.1  Synchronous  Averaging 

The  synchronous  averaging  process  for  a  continuous  vibration  signal  is  illustrated  in 
Figure  1.  Note  that  the  tachometer  must  be  of  type  that  gives  a  reference  phase  angle 
(e.g.  a  pulse  at  a  specific  shaft  angle);  a  tachometer  that  simply  gives  a  voltage 
proportional  to  speed  is  not  suitable. 


Accelerometer 


Tacho 


_n_jL 

Period 
of  shaft 
A 


^  _J  1 _ 1 1 _ 

Period  of  shaft  C 

Tacho  Multiplier 

A 
x  — 

C 

Synchronous 
average  of  shaft  C 


Figure  1.  Synchronous  Averaging. 


In  practice,  synchronous  averaging  is  usually  implemented  within  a  digital  computer 
using  discretely  sampled  vibration  and  phase  reference  signals.  This  therefore  requires 
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a  mechanism  to  ensure  the  vibration  data  are  sampled  at  equally  spaced  angular 
increments  of  the  shaft  of  interest.  That  is,  if  P  samples  are  required  for  each  revolution 
of  the  shaft  of  interest,  then  the  sample  rate  in  the  angle-domain  must  be  P/27t  for  that 
shaft,  regardless  of  any  perturbations  in  shaft  speed.  The  technique  used  at  AMRL  to 
accomplish  this  is  described  below. 

2.1.1  AMRL  Synchronous  Averaging  Technique 

1.  The  vibration  and  reference  phase  signals  are  simultaneously  sampled  at  a  rate  that 
is  sufficient  to  avoid  abasing,  and  provide  at  least  as  many  samples  per  revolution 
of  the  shaft  of  interest  as  required. 

2.  The  equivalent  cumulative  phase  angles  of  the  data  are  calculated  by: 

a)  searching  the  phase  reference  data  for  the  zero-crossings  (these  are  linearly 
interpolated  between  the  points  on  either  side  of  the  zero-crossing  level2);  and 

b)  linearly  interpolating  the  phase  angle  of  each  sample  between  the  zero-crossings. 

3.  The  vibration  data  are  then  digitally  re-sampled  with  fifth  order  splines,  using  the 
cumulative  phase  angle  results  to  provide  the  re-sampling  positions.  Fifth  order 
splines  have  been  shown  to  be  as  good  as  ideal  low-pass  filters  when  reconstructing 
a  signal  from  data  with  a  dynamic  range  less  than  75  dB  (eg  from  a  12-bit  A-D 
converter)  [2], 

4.  The  re-sampled  data  are  then  split  into  segments  corresponding  to  one  rotation 
period  of  the  shaft  of  interest,  and  ensemble-averaged. 

2.2  Effect  of  Reference  Phase  Jitter 

The  synchronous  averaging  process  can  be  modelled  in  the  angle-domain  by  the 
following  equation: 

■j  AM 

3c(#)  =  —  ^  x{6  +  nlit)  1 

™  «=0 

where: 

x  is  the  vibration  signal, 

6  is  the  shaft  angle  of  the  averaged  vibration  signal  over  the  period  [0,2;r],  and 
N  is  the  number  of  shaft  revolutions  averaged. 

If  for  each  shaft  angle,  0,  the  reference  phase  jitter  angle,  (j>,  is  modelled  as  a  random 
variable,  then  Equation  1  can  be  expressed  as: 

i  AM 

x(0)  =  —  ^x(0  +  nln  +  <f>9  (n))  2 


2  The  zero-crossing  level  need  not  necessarily  be  zero  Volts.  In  the  data  examined  here,  a  level 
mid  way  between  the  minimum  and  maximum  values  of  the  phase  reference  signal  was  used. 
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It  is  intuitive  that,  for  a  particular  level  of  phase  jitter,  the  higher  frequencies  of  the 
synchronous  average  will  be  attenuated  more  than  the  lower  frequencies.  This 
relationship  can  be  determined  by  taking  the  Fourier  Transform  of  Equation  2. 


where: 

X  is  the  Fourier  Transform  of  x,  and 
y  is  the  frequency  in  shaft  orders  (cycles  /  shaft  cycle). 


If  the  phase  jitter  variable  is  now  assumed  to  be  uniform  across  all  shaft  angles  (ie  is 
independent  of  9),  Equation  3  becomes: 

X(y)  =  1 X  ffx(/)e^+^  ejV~r)0  df  d6 

N  n-0 
™  n= 0 

_  Z(y)— 

N  n= 0 


where: 


8  is  the  Dirac  delta  function. 


It  can  be  seen  from  Equation  4  that  the  summation  is  the  combined  effect  of 
synchronous  averaging  and  phase  jitter  on  the  Fourier  Transform  of  the  vibration.  At 
the  non-synchronous  frequencies  the  summation  will  approach  zero  (ie  where  y  * 
integer),  while  at  the  synchronous  frequencies: 
ejyn2n  _  i  ^  an  integer 

The  attenuation  factor.  A,  of  the  synchronous  frequency  components  is  therefore: 
i  N-l  ,  . 

A(y)  =  —  £  ejr  ^n’  ,  y  an  integer  5 

™  „=o 

By  dropping  the  requirement  for  y  to  be  an  integer.  Equation  5  will  instead  represent 
the  envelope  of  the  attenuation  factor.  In  the  limit  asN->  <x>,  this  becomes: 

A(y)=  \ff\eiy^  d(f>  6 

where  P[$  is  the  probability  density  function  of  <f>. 
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If  the  jitter  is  assumed  to  have  a  normal  distribution  with  zero  mean  and  standard 
deviation  a,  Equation  6  becomes: 


Since,  in  practice,  N  will  be  finite,  the  actual  attenuation  factor  can  be  expected  to  lie  in 
some  band  about  Equation  7.  Since  the  difference  between  the  actual  attenuation  factor 
and  Equation  7  can  be  interpreted  as  a  variation  in  the  standard  deviation  of  the  jitter, 
the  statistical  confidence  interval  for  cr  can  be  used  to  place  confidence  limits  on  this 
band.  The  confidence  interval  for  cr[3]  is  given  by 

,  cr 

where: 

N  number  of  revolutions  averaged, 

lOO(l-a)  is  the  percentage  confidence  interval,  and 

denotes  the  chi-square  distribution  with  N- 1  degrees  of  freedom. 


For  example,  if  N  =  101,  the  99%  confidence  interval  for  the  standard  deviation,  <y, 
would  be 


1.2187c-] 


and,  as  N  increases,  the  interval  will  become  smaller.  This  is  illustrated  in  Figure  2, 
where  Equation  7  is  plotted  for  various  a,  and  the  95%  confidence  bands  for  N  - 101 
and  N  =  3993  are  shaded  about  each  curve. 


Note,  however,  that  without  an  absolute  phase  reference,  as  is  the  case  in  this 
investigation  (since  both  signals  can  be  expected  to  have  some  inherent  level  of  phase 
jitter).  Equation  7  cannot  be  used  to  determine  the  absolute  attenuation  of  each 
synchronous  average.  The  ratio  of  attenuation  factors,  however,  can  be  expressed  as: 


3  This  is  the  value  of  N  used  for  the  synchronous  averages  of  the  tail  rotor  gearbox  output  shaft 
discussed  in  the  later  sections. 
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A M  e>f 

A2  M  8 

~r  Y2  (°f  ~a2 ) 

~ei 

That  is,  the  ratio  of  the  attenuation  factors  has  the  same  form  as  the  individual  factors 
(as  long  as  a\  >  05),  with  the  variance  replaced  by  the  difference  in  variance. 

Arplitude  Response  (95%Corf.,N=  101)  Arpfrtude  Response  (95%Oxf„  N  =  399) 


Figure  2.  Attenuation  Factor  far  Various  Levels  of  Normally  Distributed  Reference  Phase  Jitter 
with  Shaded  95%  Confidence  Bands,  far  N=101  and  N=399. 


3.  AMRL  Transmission  Vibration  Monitoring  System 

3.1  S-70A-9  Transmission 

The  S-70A-9  transmission  is  shown  in  Figure  3.  The  individual  shaft  speeds  and  gear 
mesh  frequencies  are  listed  in  Appendix  1.  The  transmission  consists  of: 

a)  two  input  modules  that  are  driven  by  the  engines  at  power  turbine  speed 
(~20,900  rpm); 

b)  two  accessory  modules  that  drive  the  hydraulic  pumps  and  electrical 
generators; 

c)  a  main  module  that  drives  the  main  rotor  and  the  tail  drive  shaft; 

d)  an  intermediate  gearbox  at  the  base  of  the  tail  pylon;  and 

e)  a  tail  rotor  gearbox  at  the  top  of  the  tail  pylon. 
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Starboard 


Figure  3.  S-70A-9  Transmission. 

3.2  System  Hardware 


The  AMRL-developed  transmission  vibration  monitoring  system  is  shown 
schematically  in  Figure  4.  It  consists  of: 

a)  a  portable  ruggedised  computer  with  internally  mounted  signal  conditioning, 
anti-alias  filtering,  and  analogue-to-digital  converter  cards; 

b)  a  connector  interface  on  the  side  of  the  computer; 

c)  accelerometers  mounted  on  the  transmission  modules/  gearboxes; 

d)  a  tail  rotor  photocell  and  photocell  processor;  and 

e)  cables  to  connect  the  computer  to  the  transducers  and  aircraft  power. 
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3.3  Phase  Reference  Signals 

3.3.1  Main  Transmission 

Since  the  AC  generators  are  directly  driven  by  the  accessory  gearbox  modules  (see 
Figure  3),  a  single  phase  of  the  115  VAC  power  supply  can  be  used  as  a  phase  reference 
signal  for  the  main  transmission.  Due  to  the  construction  of  the  generator,  the  supply 
frequency  is  exactly  twice  its  shaft  speed.  Thus,  although  the  transmission  is  kept  at  a 
'constant'  speed  by  the  engine  controllers,  any  perturbation  in  speed  will  have  a 
corresponding  perturbation  in  the  supply  frequency.  From  the  data  examined  in  this 
report,  the  average  supply  frequency  was  calculated  to  be  approximately  395  Hz. 

The  signal  from  the  AC  generator  was  obtained  from  a  power  receptacle  in  the  cabin 
roof.  This  receptacle  normally  receives  power  from  the  starboard  generator.  The 
amplitude  of  the  signal  was  attenuated  from  115  VAC  to  approximately  4  VAC  by  a 
resistor  network  on  the  signal  conditioning  card. 

3.3.2  Tail  Rotor 

The  tail  rotor  phase  reference  signal  was  obtained  from  a  photocell  and  photocell 
processor  (Chadwick-Helmuth  Parts  10200  &  10170)  obtained  from  a  rotor  track  and 
balance  kit.  The  photocell  was  mounted  on  top  of  the  tail  pylon  approximately  480  mm 
behind  the  centre-line  of  the  tail  rotor  shaft.  It  was  pointed  directly  at  the  tail  rotor  disc 
with  a  slight  declination  (approx.  15°  below  the  horizontal)  to  reduce  interference  from 
direct  sunlight.  A  piece  of  reflective  tape  approximately  20  mm  x  50  mm  was  fixed  to 
the  inner  leading  edge  of  one  of  the  tail  rotor  blades.  The  beam  length  was 
approximately  275  mm. 

The  10  V  sharp  edged  pulses  from  the  photocell  processor  were  reduced  to 
approximately  1.7  V  by  a  resistor  network  on  the  signal  conditioning  card.  The  pulse 
frequency  was  approximately  20  Hz. 


4.  Flight  Data 


4.1  Data  Acquisition 

Data  from  the  two  phase  reference  signals  and  the  tail  rotor  gearbox  vibration  signal 
were  simultaneously  acquired  using  the  parameters  listed  in  Table  1  for  each  of  the 
flight  conditions  listed  in  Table  2. 


8 


DSTO-TR-0739 


Table  1.  Data  Acquisition  Parameters. 


Channels 

Tail  rotor  gearbox  accelerometer. 

Photocell  signal,  and 

AC  Generator  signal 

Total  Sample  Rate 

42,105.26  Hz  (4.00  MHz  ±0.01%  clock  /  95) 

Sample  Rate  /  Channel 

14,035.09  Hz 

Acquisition  Period 

22s 

Anti-alias  Filter 

Tail  rotor  accelerometer 

5,000  Hz 

Photocell  signal 

Bypass 

AC  Generator  signal 

Bypass 

Table  2.  Flight  Conditions. 


# 

Flight  Condition 

Altitude  (ft) 

KIAS 

Engine  Torque  (% /  %) 

Nr  (%) 

A 

Straight  &  Level  60% 

1300 

110 

60/60 

100 

B 

Straight  &  Level  70% 

1300 

120 

70/70 

100 

C 

Straight  &  Level  80% 

1100 

125 

80/80 

100 

D 

Straight  &  Level  90% 

1100 

135 

90/90 

100 

E 

Straight  &  Level  100% 

1100 

140 

100/100 

100 

F 

Right  Turn  -  30°  Bank 

1300 

120 

80/80 

100 

G 

Right  Turn  -  45°  Bank 

1300 

120 

80/80 

100 

H 

Hover 

25 

0 

60/60 

100 

I 

Clockwise  Tail  Rotor  Turn 

25 

0 

60/60 

100 

J 

Anti-clockwise  Tail  Rotor  Turn 

25 

0 

60/60 

100 

4.2  Data  Processing 

The  raw  data  files  were  first  processed  with  AMRL  software,  and  the  results  were  then 
analysed  using  MATLAB  (see  Appendix  2  for  the  M  files),  as  described  below. 

4.2.1  Raw  Data 

For  each  flight  condition: 

1.  The  equivalent  cumulative  phase  angles  of  the  data  were  calculated  (see  §2.1.1). 

2.  A  separate  synchronous  average  (see  Table  3)  of  the  tail  rotor  gearbox  output 
shaft  vibration  was  calculated  using  the  each  of  the  two  different  phase  reference 
data  records. 


Table  3.  Synchronous  Average  Parameters  for  Tail  Rotor  Gearbox  Output  Shaft 


Re-sample  Rate 
(Samples/ Rev) 

Shaft  Revs 
Averaged 

Ratio  of  Shaft  to  AC  Generator  Signal 

Ratio  of  Shaft  to 
Photocell  Signal 

1024 

3994 

1/2  X  37/76  X  21/100  *  75/22  x  25/31  x  19/53 
=  1107225  /  21976768 

i/i 

4  Chosen  as  a  multiple  of  the  number  of  teeth  on  the  mating  pinion  (19)  to  minimise  the 
modulation  of  the  output  gear  average  by  the  input  pinion. 
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4.2.2  Rotor  Speed  Perturbations 

While  the  main  rotor  speed  in  the  Black  Hawk  is  nominally  constant  (the  engines  are 
controlled  to  a  constant  speed),  it  is  inevitable  that  load  fluctuations  will  induce  small 
perturbations  about  this  speed.  Since  a  constant  speed  implies  that  the  phase  angle 
should  increase  linearly,  these  perturbations  can  be  examined  by  subtracting  a  straight 
line  fit  from  the  cumulative  phase  data.  To  enable  the  same  fit  to  be  used  for  both  data 
series,  the  AC  generator  data  were  first  converted  to  tail  rotor  phase  angle  units  by: 

a)  multiplying  the  AC  generator  data  by  the  gear  ratio  between  the  generator  and 
the  tail  rotor  (1107225/21976768), 

b)  performing  a  linear  regression  on  the  data  using  the  MATLAB  polyfit  and 
polyval  functions, 

c)  subtracting  this  linear  fit  from  both  signals,  and 

d)  removing  any  offset  between  the  data  series  by  subtracting  their  means. 

4.2.3  Relative  Phase  Difference 

The  relative  difference  (jitter)  between  the  cumulative  phase  data  were  examined  by 
converting  the  AC  generator  data  into  tail  rotor  phase  angle  units,  then  subtracting 
these  data  from  the  photocell  data,  and  estimating  the  probability  density  function  of 
the  difference.  This  was  done  by: 

a)  multiplying  the  AC  generator  data  by  the  gear  ratio  between  the  generator  and 
the  tail  rotor  (1107225/21976768), 

b)  subtracting  the  AC  generator  data  from  the  photocell  data, 

c)  removing  any  zero-offset  of  the  difference  data  by  subtracting  the  mean, 

d)  converting  the  difference  data  to  degrees,  and 

e)  calculating  a  histogram  of  the  difference  data  using  the  MATLAB  hist  function. 

4.2.4  Synchronous  Average  Difference 

To  compare  the  synchronous  averages,  they  were  first  aligned  to  remove  any  phase 
offset5,  then  the  AC  generator  average  was  subtracted  from  the  photocell  average.  This 
was  done  by: 

a)  finding  the  phase  angle  corresponding  to  the  point  of  maximum  (circular) 
cross-correlation  between  the  averages,  using  the  MATLAB  fmin  function, 

b)  phase-shifting  the  AC  generator  average  by  this  angle,  and 

c)  subtracting  the  phase-shifted  AC  generator  average  from  the  photocell  average. 


5  The  phase  offsets  between  the  synchronous  averages  were  an  artefact  of  the  averaging  process, 
which  was  arbitrarily  started  at  the  first  zero-crossing  of  the  respective  phase  reference 
signals. 


10 


DSTO-TR-0739 


5.  Results 

The  results  from  the  rotor  speed  perturbations,  relative  phase  difference,  probability 
density  function  estimates,  and  synchronous  average  difference  analyses  are  presented 
below. 

5.1  Rotor  Speed  Perturbations 

The  effects  of  rotor  speed  perturbations  on  the  tail  rotor  phase  angle  are  shown  from 
Figure  5  to  Figure  14  (see  §4.2.2  for  description  of  method).  The  vertical  axis  shows  the 
deviation  of  the  tail  rotor  angle  from  that  which  would  be  expected  if  the  tail  rotor 
speed  was  exactly  constant.  In  each  figure,  the  plot  on  the  left  shows  the  entire  data 
record,  while  the  plot  on  the  right  shows  a  selected  segment. 

The  following  observations  can  be  made  about  these  figures: 

1.  The  rotor  speed  perturbations  during  forward  flight  (straight  and  level,  and  banked 
turns)  are  smaller  than  those  during  hover  flight  (hover,  and  tail  rotor  turns).  This  is 
evidence  of  larger  torque  variations  during  hover,  probably  caused  by  the  pilot 
trying  to  maintain  station  at  a  fixed  point  over  the  ground. 

2.  On  the  macro  scale,  the  left  hand  plots  reveal  virtually  no  detectable  separation 
between  the  photocell  and  AC  generator  phase  data,  showing  that  both  signals  seem 
to  follow  each  other  very  well,  which  is  as  would  be  expected  for  a  geared 
drivetrain. 

3.  On  the  micro  scale,  however,  the  right  hand  plots  reveal  that  there  are  detectable 
differences  between  the  phase  data  which  can  be  characterised  by  their  flight 
conditions. 

a)  Forward  flight:  The  photocell  data  appears  to  have  a  small  degree  of  "overshoot" 
compared  to  the  AC  generator  data.  It  can  be  seen  that  where  the  curves  change 
direction,  corresponding  to  an  acceleration  or  deceleration  of  the  rotors,  the  AC 
generator  data  responds  before  tire  photocell  data.  This  is  probably  because  the 
phase  angles  are  calculated  with  a  zero-crossing  algorithm  (see  §2.1.1),  and  the 
AC  generator  signal  has  a  higher  frequency  than  the  photocell  signal  (395  Hz 
compared  to  20  FIz). 

b)  Hover  Flight:  In  these  cases,  there  are  larger  and  longer  separations  between  the 
photocell  and  AC  generator  data.  While  this  may  still  be  partially  attributable  to 
the  "overshoot"  mentioned  above,  the  sustained  separation  of  the  data  is  more 
likely  to  be  due  to  twisting  of  the  long  tail  drive  shaft  caused  by  variation  of  the 
tail  rotor  torque.  The  presence  and  the  characteristic  duration  of  these  deviations 
is  consistent  with  pedal  control  movements  made  by  the  pilot,  which  would  be 
expected  to  be  greater  during  hover  and  tail  rotor  turns,  than  in  forward  flight. 
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Tail  Rotor  Angle:  Deviation  from  Constant  Speed  Value  Tail  Rotor  Angle:  Deviation  from  Constant  Speed  Value 


Time  (s)  Time  (s) 

.  Effect  of  Rotor  Speed  Perturbations  on  the  Tail  Rotor  Phase  Angle  for  Straight  and 
Level  Flight,  60%  Torque. 

Tail  Rotor  Angle:  Deviation  from  Constant  Speed  Value  Tail  Rotor  Angle:  Deviation  from  Constant  Speed  Value 


Time  (s)  Time  (s) 

Figure  6.  Effect  of  Rotor  Speed  Perturbations  on  the  Tail  Rotor  Phase  Angle  for  Straight  and 
Level  Flight,  70%  Torque. 
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Tail  Rotor  Angle:  Deviation  from  Constant  Speed  Value  Tail  Rotor  Angle:  Deviation  from  Constant  Speed  Value 


Time  (s)  Time  (s) 

Effect  of  Rotor  Speed  Perturbations  on  the  Tail  Rotor  Phase  Angle  for  Straight  and 
Level  Flight,  80%  Torque. 

Tail  Rotor  Angle:  Deviation  from  Constant  Speed  Value  Tail  Rotor  Angle:  Deviation  from  Constant  Speed  Value 


Time  (s)  Time  (s) 

Figure  8.  Effect  of  Rotor  Speed  Perturbations  on  the  Tail  Rotor  Phase  Angle  for  Straight  and 
Level  Flight,  90%  Torque. 
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Tail  Rotor  Angle:  Deviation  from  Constant  Speed  Value  Tail  Rotor  Angle:  Deviation  from  Constant  Speed  Value 


Time  (s)  Time  (s) 

Effect  of  Rotor  Speed  Perturbations  on  the  Tail  Rotor  Phase  Angle  for  Straight  and 
Level  Flight,  100%  Torque. 

Tail  Rotor  Angle:  Deviation  from  Constant  Speed  Value  Tail  Rotor  Angle:  Deviation  from  Constant  Speed  Value 


Time  (s)  Time  (s) 

Figure  10.  Effect  of  Rotor  Speed  Perturbations  on  the  Tail  Rotor  Phase  Angle  for  Right  Turn 
30°  Bank. 
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Tail  Rotor  Angle:  Deviation  from  Constant  Speed  Value  Tail  Rotor  Angle:  Deviation  from  Constant  Speed  Value 


Figure  11.  Effect  of  Rotor  Speed  Perturbations  on  the  Tail  Rotor  Phase  Angle  for  Right  Turn , 
45°  Bank. 

Tail  Rotor  Angle:  Deviation  from  Constant  Speed  Value  Tail  Rotor  Angle:  Deviation  from  Constant  Speed  Value 


Figure  12.  Effect  of  Rotor  Speed  Perturbations  on  the  Tail  Rotor  Phase  Angle  for  Hover. 
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Tail  Rotor  Angle:  Deviation  from  Constant  Speed  Value  Tail  Rotor  Angle:  Deviation  from  Constant  Speed  Value 


Figure  13.  Effect  of  Rotor  Speed  Perturbations  on  the  Tail  Rotor  Phase  Angle  for  Clockwise  Tail 
Rotor  Turn. 


Tail  Rotor  Angle:  Deviation  from  Constant  Speed  Value 


Tail  Rotor  Angle:  Deviation  from  Constant  Speed  Value 


Figure  14.  Effect  of  Rotor  Speed  Perturbations  on  the  Tail  Rotor  Phase  Angle  for  Anti¬ 
clockwise  Tail  Rotor  Turn. 

5.2  Relative  Phase  Difference 

The  relative  differences  between  the  two  tail  rotor  phase  data  sets  for  each  flight 
condition  (obtained  by  subtracting  one  data  set  from  the  other)  are  shown  from  Figure 
15  to  Figure  24.  The  most  distinguishing  feature  of  these  plots  is  the  disparity  between 
results  for  the  forward  flight  and  hover  flight  conditions.  Apart  from  the  seemingly 
random  high  frequency  'noise'  present  in  the  forward  flight  results,  these  plots  are 
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relatively  flat,  with  only  a  few  underlying  deviations  evident  beneath  the  noise.  On  the 
other  hand,  the  hover  flight  results  exhibit  quite  large  underlying  deviations. 
Generally,  these  features  correlate  well  with  the  phase  separations  observed  in  §5.1. 
From  this,  it  is  postulated  that,  in  general: 

a)  the  underlying  large  scale  deviations  correspond  to  the  dynamic  twisting  of  the 
tail  drive  shaft  brought  about  by  pilot  induced  tail  rotor  torque  variations,  and 

b)  the  'noise'  is  most  probably  attributable  to  measurement  error  (see  §5.5). 


Phase  Difference 


Figure  15.  Relative  Phase  Difference  for 
Straight  &  Level  Flight,  60% 
Torque. 


Phase  Difference 


Figure  16.  Relative  Phase  Difference  for 
Straight  &  Level  Flight,  70% 
Torque. 


Phase  Difference 


Figure  17.  Relative  Phase  Difference  for 
Straight  &  Level  Flight,  80% 
Torque. 


Phase  Difference 


Figure  18.  Relative  Phase  Difference  for 
Straight  &  Level  Flight,  90% 
Torque. 
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Phase  Difference 


Figure  19.  Relative  Phase  Difference  for 
Straight  &  Level  Flight,  100% 
Torque. 


Phase  Difference 


Figure  20.  Relative  Phase  Difference  for  Right 
Turn,  SO0  Bank. 
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Figure  21.  Relative  Phase  Difference  for  Right 
Turn,  45°  Bank. 


Phase  Difference 


Figure  22.  Relative  Phase  Difference  for 
Hover. 


Phase  Difference 


Figure  23.  Relative  Phase  Difference  for 
Clockwise  Tail  Rotor  Turn. 


Phase  Difference 


Figure  24.  Relative  Phase  Difference  for  Anti¬ 
clockwise  Tail  Rotor  Turn. 
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5.3  Probability  Density  Function  Estimates 

The  probability  density  function  (PDF)  estimates  of  the  relative  phase  difference  data 
from  §5.2  are  shown  from  Figure  25  to  Figure  34.  They  were  calculated  across  a  band 
from  -1.33°  to  1.33°,  with  133  bins  of  width  0.020°.  The  standard  deviation,  minimum 
and  maximum  are  shown  in  the  upper  right  comer  of  each  plot.  For  comparison,  a 
normal  distribution  with  the  same  standard  deviation  has  been  overlaid  on  top  of  each 
PDF. 

It  can  be  seen  from  these  plots  that  the  distinguishing  features  previously  identified  in 
the  phase  difference  data  (see  §5.2)  are  directly  translated  into  their  probability  density 
functions.  The  plots  can  therefore  be  split  into  the  same  two  groups:  forward  flight,  and 
hover  flight.  The  following  observations  can  be  made  about  these  two  groups. 

5.3.1  Forward  Flight 

1.  Each  PDF  estimate  very  closely  approximates  the  overlaid  normal  distribution.  The 
small  departures  from  the  normal  distribution  curve  can  be  attributed  to  the  small 
underlying  deviations  evident  in  the  corresponding  figures  of  §5.2.  Note,  however, 
that  even  though  the  difference  data  approach  a  normal  distribution,  this  does  not 
necessarily  imply  that  the  underlying  jitter  inherent  in  each  phase  reference  signal  is 
also  normally  distributed,  although  it  is  likely. 

2.  The  standard  deviations  of  the  distributions  decrease  slightly  with  engine  torque, 
from  0.139°  at  60%  torque,  to  0.102°  at  100%  torque.  This  would  appear  to  correlate 
well  with  the  magnitude  of  the  underlying  deviations  observed  in  the 
corresponding  figures  of  §5.2.  From  this,  and  the  postulation  that  they  were  caused 
by  tail  rotor  torque  variation,  it  can  be  inferred  that  this  decreases  slightly  as  the 
forward  airspeed  increases  (from  110  to  140  KIAS). 

5.3.2  Hover  Flight 

1.  These  distributions  do  not  match  the  overlaid  normal  distributions.  The  large 
departures  from  the  normal  curve  can  be  attributed  to  the  large  underlying 
deviations  caused  by  the  dynamic  tail  drive  shaft  twist  evident  in  the  corresponding 
figures  of  §5.2. 

2.  The  standard  deviations  of  the  distributions  increase  for  the  tail  rotor  turns,  from 
0.204°  for  hover,  to  0.441°  for  the  clockwise  turn,  and  0.383°  for  the  anticlockwise 
turn.  This  is  probably  due  to  the  pilot  making  larger  pedal  adjustments  during  the 
tail  rotor  turns,  than  while  hovering. 
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PDF  of  Phase  Difference  PDF  of  Phase  Difference 


Figure  25.  PDF  Estimate  of  Phase  Difference 
for  Straight  &  Level  Flight,  60% 
Torque. 


Figure  28.  PDF  Estimate  of  Phase  Difference 
for  Straight  &  Level  Flight,  90% 
Torque. 


Figure  26.  PDF  Estimate  of  Phase  Difference 
for  Straight  &  Level  Right,  70% 
Torque. 


Figure  29.  PDF  Estimate  of  Phase  Difference 
for  Straight  &  Level  Flight,  100%> 
Torque. 


PDF  of  Phase  Difference  PDF  of  Phase  Difference 


Figure  27.  PDF  Estimate  of  Phase  Difference  Figure  30.  PDF  Estimate  of  Phase  Difference 
for  Straight  &  Level  Flight,  80%  for  Right  Turn,  30°  Bank. 

Torque. 
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PDF  of  Phase  Difference  PDF  of  Phase  Difference 


Figure  31.  PDF  Estimate  of  Phase  Difference 
for  Right  Turn,  45°  Bank. 


Figure  32.  PDF  Estimate  of  Phase  Difference 
for  Hover. 


Figure  33.  PDF  Estimate  of  Phase  Difference 
for  Clockwise  Tail  Rotor  Turn. 


PDF  of  Phase  Difference 


Figure  34.  PDF  Estimate  of  Phase  Difference 
for  Anti-Clockwise  Tail  Rotor 
Turn. 


5.4  Synchronous  Averages 

The  synchronous  averages  of  the  tail  rotor  gearbox  output  shaft  vibration  are  shown 
from  Figure  39  to  Figure  48.  Each  figure  is  split  into  three  sections  showing: 

a)  the  synchronous  average  and  its  spectrum  computed  using  the  photocell  signal, 

b)  the  synchronous  average  and  its  spectrum  computed  using  the  AC  generator 
signal,  and 

c)  the  difference  between  the  synchronous  averages,  and  the  amplitude  ratios  of 
the  gear  mesh  harmonics. 

For  comparison  purposes,  the  signal  amplitudes  have  been  extracted  from  these  figures 
and  summarised  in  Table  4. 

Before  discussing  these  results,  however,  the  following  information  in  §5.4.1  and  §5.4.2 
should  be  noted. 
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5.4.1  Attenuation  due  to  Reference  Phase  Jitter 

The  theoretical  attenuation  of  a  synchronous  average  due  to  reference  phase  jitter  was 
developed  in  §2.2.  Application  of  this  theory  to  die  data,  however,  is  limited  by  two 
factors: 

1.  Firstly,  there  is  no  absolute  measure  of  the  phase  jitter  in  each  phase  reference 
signal,  as  both  must  be  considered  to  have  at  least  some  measurement  error  (see 
§5.5).  Also,  the  tail  rotor  photocell  signal  cannot  automatically  be  assumed  to 
provide  a  more  accurate  measure  of  die  tail  rotor  angle  than  the  AC  generator 
signal.  In  fact,  the  higher  frequency  of  the  AC  generator  signal  (395  Hz),  compared 
to  the  photocell  signal  (20  Hz),  may  allow  it  to  better  follow  any  short  term  speed 
perturbations  (disregarding  dynamic  shaft  twisting)  (see  §5.1). 

Without  an  absolute  phase  reference,  the  absolute  attenuation  of  the  synchronous 
vibration  frequencies  cannot  be  determined.  All  that  can  be  examined  is  the  ratio  of 
the  attenuation  factors.  From  Equation  8,  this  was  shown  to  have  the  same  form  as 
the  individual  attenuation  factors  for  each  average  (as  long  as  oi  >  05),  but  with  the 
jitter  variance  replaced  by  the  difference  in  jitter  variance.  If  the  assumptions  made 
in  deriving  Equation  8  hold,  then  the  amplitude  ratios  of  the  gear  mesh  harmonics 
should  follow  the  general  shape  of  the  curves  shown  in  Figure  2. 

2.  Secondly,  as  can  be  seen  in  §5.3,  the  assumption  in  §2.2  of  normally  distributed 
phase  jitter  is  not  necessarily  valid  for  any  of  the  data  examined  here.  It  would 
appear  that  the  measurement  error  component  of  the  jitter  comes  close  to  being 
normally  distributed,  but  the  tail  drive  shaft  twist  component  does  not.  Attempts 
were  made  to  separate  the  measurement  error  from  the  underlying  dynamic  tail 
drive  shaft  twist  (so  that  they  might  be  investigated  separately)  by  filtering  the  data, 
but  this  was  unsuccessful  because  no  matter  where  the  cross-over  frequency  was 
placed,  some  of  the  measurement  error  would  be  left  mixed  in  with  the  shaft  twist, 
and  vice  versa. 

Since  the  assumption  of  normally  distributed  phase  jitter  does  not  necessarily  apply, 
the  amplitude  ratios  of  the  gear  mesh  harmonics  cannot  be  expected  to  follow 
exactly  the  curves  shown  in  Figure  2.  The  more  the  jitter  departs  from  a  normal 
distribution,  the  more  the  amplitude  ratios  can  be  expected  to  depart  from  these 
curves. 

5.4.2  Extraneous  Vibration  Impulses 

Large  impulses  were  observed  in  the  tail  rotor  gearbox  vibration  data  for  all  the 
forward  flight  conditions,  but  not  during  the  hover  flight  conditions.  An  example  of 
these  impulses  is  shown  in  Figure  35,  with  the  detail  of  a  typical  impulse  shown  in 
Figure  36.  Although  impulses  are  typically  symptomatic  of  a  loose  accelerometer,  this 
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was  found  not  to  be  the  case.  It  is  therefore  speculated  that,  since  the  impulses  only 
occurred  during  forward  flight,  they  were  probably  caused  by  the  aerodynamic 
interaction  between  the  tail  rotor  blades  and  the  main  rotor  tip  vortices.  That  is,  during 
forward  flight,  the  main  rotor  tip  vortices  were  being  swept  back  through  the  tail  rotor 
disc,  and  wherever  they  intersected  with  a  tail  rotor  blade  they  would  generate  a  shock 
in  the  blade  that  would  travel  through  the  hub  to  the  gearbox.  During  hover  the 
vortices  would  descend  vertically  away  from  the  main  rotor  and  not  interact  with  the 
tail  rotor.  If  this  speculation  is  correct,  then  the  impulses  are  not  desirable  for  health 
monitoring  because  they  may  mask  faults  within  the  gearbox. 

It  is  difficult  to  determine  the  extent  to  which  the  impulses  have  affected  the 
synchronous  averages  shown  here.  As  the  impulses  were  slightly  non-synchronous, 
they  have  been  attenuated  by  the  averaging  process,  leaving  little  visual  evidence  of 
their  presence  in  the  averages.  This  is  not  to  say  that  their  effects  have  been  entirely 
removed,  however,  as  their  presence  alone  has  forced  the  expansion  of  the  ADC  input 
range,  effectively  reducing  dynamic  range  of  the  gear  mesh  vibration  data.  For  the 
purposes  of  this  investigation,  however,  any  remaining  effects  should  not  affect 
comparisons  between  the  averages  greatly,  since  the  same  vibration  data  were  used  for 
each  average. 


Figure  35.  Typical  Tail  Rotor  Gearbox  Figure  36.  Typical  Impulse. 
Vibration  During  Forward  Flight. 


5.4.3  Forward  Flight 

The  following  observations  can  be  made  about  the  forward  flight  synchronous 
averages. 

1.  In  general,  and  especially  when  compared  to  hover  flight  (see  §5.4.4),  there  is  good 
overall  agreement  between  the  AC  generator  and  photocell  averages.  This  indicates 
that  where  dynamic  twisting  of  the  tail  drive  shaft  is  minimised,  both  phase 
reference  signals  perform  similarly  well  for  synchronous  averaging  purposes. 
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2.  It  can  be  seen  from  Table  4,  however,  that  the  overall  amplitudes  of  the  AC 
generator  averages  are,  in  all  but  one  case  (flight  condition  E  - 100%  torque),  slightly 
larger  than  the  photocell  averages.  This  indicates,  in  a  general  sense,  that  the  AC 
generator  may  provide  a  fractionally  more  accurate  phase  reference  signal  in  these 
flight  conditions  than  the  photocell,  probably  because  of  its  higher  frequency  and 
lack  of  "overshoot"  (see  §5.1).  These  differences  are  quite  small,  however,  and  do 
not  necessarily  justify  the  preferential  selection  of  the  AC  generator  signal  over  the 
photocell  signal. 

3.  The  amplitudes  of  the  gear  mesh  harmonics  (multiples  of  53  shaft  orders)  in  the  AC 
generator  averages  are  nearly  all  slightly  greater  than,  or  equal  to,  those  in  the 
photocell  averages.  This  provides  some  additional  evidence  that  the  AC  generator 
signal  might  be  a  fractionally  better  phase  reference  than  the  photocell  signal.  It 
should  be  noted,  however,  that  the  amplitude  ratios,  plotted  in  Figure  37,  for  the 
first  four  gear  mesh  harmonics  are  all  very  close  to  a  value  of  one,  which  is  the 
expected  ratio  for  phase  reference  signals  with  equal  levels  of  jitter  (see  §5.4.1). 
These  results  therefore  also  show  some  qualitative  agreement  with  the  mathematical 
model.  The  larger  amplitude  variations  for  the  fifth  harmonic  may  be  caused  by  this 
frequency  falling  in  the  transition-band  of  the  anti-alias  filter,  as  small  changes  in 
rotor  speed  will  have  affected  the  attenuation  of  this  frequency  differently  to  the 
other  harmonics  (due  to  the  slope  of  the  filter  response). 

5.4.4  Hover  Flight 

The  following  observations  can  be  made  about  the  hover  flight  synchronous  averages. 

1.  In  general,  there  is  poor  overall  agreement  between  the  AC  generator  and  photocell 
averages.  This  indicates  that  the  dynamic  tail  drive  shaft  twist  observed  in  §5.2  has  a 
very  significant  effect  on  the  synchronous  averaging  of  the  tail  rotor  gearbox 
vibration. 

2.  It  can  be  seen  from  Table  4  that  the  overall  amplitudes  of  the  photocell  averages  are, 
in  two  out  of  three  cases  (hover  and  anti-clockwise  tail  rotor  turn),  greater  than  the 
AC  generator  averages.  This  provides  a  general  indicator  that  the  photocell  is  a 
more  accurate  phase  reference  in  these  flight  conditions  than  the  AC  generator.  The 
exception  for  the  clockwise  tail  rotor  turn  data  may  be  due  to  the  relatively  low 
levels  of  tail  rotor  torque  in  this  flight  condition. 

3.  With  the  exception  of  the  clockwise  tail  rotor  turn,  the  amplitudes  of  the  gear  mesh 
harmonics  in  the  photocell  averages  are  nearly  all  greater  than  those  in  the  AC 
generator  averages.  This  provides  additional  evidence  that  the  photocell  is  a  better 
phase  reference  than  the  AC  generator  in  these  flight  conditions.  The  amplitude 
ratios  for  the  gear  mesh  harmonics,  plotted  in  Figure  38,  show  very  little  agreement 
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with  the  mathematical  model  however.  As  predicted,  this  is  most  probably  due  to 
the  assumption  in  the  model  that  the  phase  jitter  is  normally  distributed. 


Table  4.  Synchronous  Average  Amplitudes 


# 

AC  Generator 

Average  (Grms) 

Photocell  Average 
(Grms) 

Difference  Average 
(Grms) 

Difference  /  AC 
Generator  (%) 

A 

4.125 

3.915  . 

0.278 

6.74 

B 

4.789 

4.656 

0.194 

4.05 

C 

4.691 

4.612 

0.156 

3.33 

D 

4.047 

4.005 

0.148 

3.66 

E 

2.756  | 

2.767 

0.112 

4.06 

F 

5.252 

5.127 

0.183 

3.48 

G 

5.137 

5.001 

0.227 

4.42 

H 

3.309 

3.834 

0.759 

22.9 

I 

4.485 

4.165 

0.801 

17.9 

J 

2.522 

3.650 

1.588 

63.0 

Figure  37.  Amplitude  Ratios  of  Gear  Mesh  Figure  38.  Amplitude  Ratios  of  Gear  Mesh 
Harmonics  for  Forward  Flight  Harmonics  for  Hover  Flight 

(Photocell /  AC  Generator).  (Photocell / AC  Generator). 
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Synchronous  Average  Spectrum 
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Figure  39.  Synchronous  Averages  for  Straight  &  Level  Flight,  60%  Torque 

a)  using  the  photocell  signal, 

b)  using  the  AC  generator  signal  (</>  = phase  offset  of(b)to  (a)), 

c)  difference  between  (a)  and  phase  shifted  (b),  and  amplitude  ratios  of  gear 
mesh  harmonics. 
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Figure  40.  Synchronous  Averages  for  Straight  &  Level  Flight,  70%  Torque 

a)  using  the  photocell  signal, 

b)  using  the  AC  generator  signal  (</>  =  phase  offset  of  (b)  to  (a)), 

c)  difference  between  (a)  and  phase  shifted  (b),  and  amplitude  ratios  of  gear 
mesh  harmonics. 
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Spectrum 
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Figure  41.  Synchronous  Averages  for  Straight  &  Level  Flight ,  80%  Torque 

a)  using  the  photocell  signal, 

b)  using  the  AC  generator  signal  (<f>  =  phase  offset  of(b)  to  (a)), 

c)  difference  between  (a)  and  phase  shifted  (b),  and  amplitude  ratios  of  gear 
mesh  harmonics. 
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Figure  42.  Synchronous  Averages  for  Straight  &  Level  Flight,  90%  Torque 

a)  using  the  photocell  signal, 

b)  using  the  AC  generator  signal  (<j>=  phase  offset  of  (b)  to  (a)), 

c)  difference  between  (a)  and  phase  shifted  (b),  and  amplitude  ratios  of  gear 
mesh  harmonics. 
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Figure  43.  Synchronous  Averages  for  Straight  &  Level  Flight,  100%  Torque 

a)  using  the  photocell  signal, 

b)  using  the  AC  generator  signal  (tj>  =  phase  offset  of  (b)  to  (a)), 

c)  difference  between  (a)  and  phase  shifted  (b),  and  amplitude  ratios  of  gear 
mesh  harmonics. 
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Figure  44.  Synchronous  Averages  for  Right  Turn,  30°  Bank 

a)  using  the  -photocell  signal, 

b)  using  the  AC  generator  signal  (<f>  =  phase  offset  of  (b)  to  (a)), 

c)  difference  between  (a)  and  phase  shifted  (b),  and  amplitude  ratios  of  gear 
mesh  harmonics. 
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Figure  45.  Synchronous  Averages  for  Right  Turn,  45°  Bank 

a)  using  the  photocell  signal, 

b)  using  the  AC  generator  signal  (<f>  =  phase  offset  of(b)  to  (a)), 

c)  difference  between  (a)  and  phase  shifted  (b),  and  amplitude  ratios  of  gear 
mesh  harmonics. 
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Figure  46.  Synchronous  Averages  for  Hover 

a)  using  the  photocell  signal, 

b)  using  the  AC  generator  signal  (</>  =  phase  offset  of  (b)  to  (a)), 

c)  difference  between  (a)  and  phase  shifted  (b),  and  amplitude  ratios  of  gear 
mesh  harmonics. 
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Figure  47.  Synchronous  Averages  for  Clockwise  Tail  Rotor  Turn 

a)  using  the  photocell  signal, 

b)  using  the  AC  generator  signal  (<j>  =  phase  offset  of  (b)  to  (a)), 

c)  difference  between  (a)  and  phase  shifted  (b),  and  amplitude  ratios  of  gear 
mesh  harmonics. 
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Figure  48.  Synchronous  Averages  for  Anti-Clockwise  Tail  rotor  Turn 

a)  using  the  photocell  signal, 

b)  using  the  AC  generator  signal  (<f)  =  phase  offset  of  (b)  to  (a)), 

c)  difference  between  (a)  and  phase  shifted  (b),  and  amplitude  ratios  of  gear 
mesh  harmonics. 
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5.5  Measurement  Errors 

Errors  in  the  measurement  of  the  signal  phase  can  be  attributed  to  the  causes  listed 
below.  It  is  not  possible  to  accurately  determine  the  overall  measurement  error, 
however,  as  this  would  require  a  measurement  system  that  is  more  accurate  than 
those  used.  Nevertheless,  the  trends  identified  in  the  phase  difference  distributions 
suggest  that  the  measurement  error  must  be  small  enough  to  produce  the  narrowest 
PDF  (straight  and  level  flight  at  100%  torque),  which  has  a  standard  deviation  of  only 
0.104°.  Thus,  the  measurement  error  would  be  expected  to  be  of  a  similar  magnitude, 
or  approximately  0.1°. 

5.5.1  Zero-Crossing  Error 

As  the  positions  of  the  zero-crossings  are  determined  by  linearly  interpolating 
between  adjacent  samples  above  and  below  the  zero-crossing  threshold,  there  will 
inevitably  be  some  difference  between  the  interpolated  zero-crossings  and  the  true 
zero-crossings.  This  error  will  diminish  as  the  sample  rate  increases. 

Typical  examples  of  the  AC  generator  and  photocell  signals,  are  shown  below 
(Figure  49  to  Figure  52).  In  both  cases,  it  can  be  seen  that  the  sample  rate  is 
sufficiently  high  to  reduce  this  error  to  a  very  low  level,  and  that  linear  interpolation 
is  not  an  unreasonable  method  of  determining  the  positions  of  the  zero-crossings.  In 
either  case,  the  maximum  error  will  be  limited  to  one  sample  period  (71.25  ps),  but 
should  be  considerably  less  than  this.  At  an  average  tail  rotor  speed  of  1190  rpm,  one 
sample  period  corresponds  to  approximately  0.5°. 


Figure  49.  Typical  AC  Generator  Signal. 


Figure  50.  Typical  Photocell  Signal. 


Figure  51.  AC  Generator  Signal  Detail. 


Seconds 


Figure  52.  Photocell  Signal  Detail. 
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5.5.2  Phase  Interpolation  Error  Between  Zero-Crossings 

As  the  phase  of  the  tail  rotor  is  linearly  interpolated  between  the  positive-going  zero- 
crossings  of  the  phase  signals,  there  will  inevitably  be  some  difference  between  the 
interpolated  phase  and  the  true  phase.  That  is,  the  true  phase  may  behave  in  a  non¬ 
linear  way  between  the  zero-crossings.  This  error  will  be  smaller  for  the  AC 
generator  signal,  as  it  has  approximately  19.85  zero-crossings  for  every  one  of  the 
photocell  signal.  The  error  for  the  photocell  could  possibly  be  reduced,  by  placing 
reflective  tape  on  the  other  tail  rotor  blades  to  give  more  pulses  per  revolution. 
However,  this  would  require  an  evaluation  of  the  lead/ lag  error  between  the  blades. 

Unfortunately,  it  is  not  possible  to  estimate  this  error  without  a  phase  reference 
signal  that  is  more  accurate  than  those  available  (eg  from  a  shaft  encoder).  The 
substantial  inertia  of  the  tail  rotor,  however,  will  act  to  limit  the  rate  of  change  of 
phase,  and  in  the  absence  of  large  impulsive  torques,  the  magnitude  of  this  error 
would  be  expected  to  be  similar  to  the  zero-crossing  error. 

5.5.3  Other  Errors 

The  following  errors  will  all  contribute  to  the  overall  measurement  error,  but  were 
considered  to  be  less  significant  than  those  above. 

1.  Photocell  Triggering:  The  exact  point  at  which  the  photocell  is  triggered  by  the 
reflective  tape  (ie  how  far  the  reflective  tape  penetrates  the  beam  before  triggering 
the  photocell)  may  vary  slightly  from  one  revolution  of  the  tail  rotor  to  the  next. 
This  will  be  compounded  to  some  extent  by  flapping  of  the  rotor  blade  under  the 
control  inputs  and  dynamic  loads  experienced  in  flight. 

2.  Signal  Noise/Interference:  Noise  or  interference  has  the  potential  to  make  the 
signal  zero-crossings  appear  earlier  or  later  than  their  true  positions.  Noise 
introduced  through  tape  recording  and  filtering  processes  has  been  found  to  be 
significant  in  previous  investigations  [4],  but  has  been  kept  to  a  minimum  here 
through  the  direct  digital  acquisition  of  data  using  the  computer.  Interference  by 
the  tail  pylon  strobes  was  observed  in  the  photocell  signal,  but  these  strobes  were 
turned  off  during  the  acquisition  of  all  data. 

3.  Analogue-to-Digital  Converter  (ADQ  Error:  ADC  errors  will  also  make  the  signal 
zero-crossings  appear  earlier  or  later  than  their  true  positions.  A  12-bit  ADC  will 
have  at  least  ±Vz  bit  of  uncertainty  in  the  measurement.  With  an  input  range  of  ±5 
V,  this  corresponds  to  an  uncertainty  of  approximately  1  mV. 

4.  AC  Signal  Distortion:  The  AC  generator  supplies  power  to  a  number  of  devices  in 
the  aircraft  all  of  which  have  the  potential  to  disturb  or  distort  the  signal 
waveform.  Localised  disturbances  that  cause  the  zero-crossing  to  move  relative  to 
its  undisturbed  position  will  affect  the  results.  Continuous  distortion  of  the 
waveform,  however,  should  not  affect  the  results,  as  the  zero-crossing  error  will 
be  constant. 
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6.  Conclusions 


This  report  has  examined  the  relative  jitter,  due  to  dynamic  tail  drive  shaft  twist, 
between  phase  reference  signals  obtained  from  the  main  and  tail  rotor  gearboxes  in  a 
S-70A-9  Black  Hawk  helicopter.  The  results  from  this,  and  a  comparison  of  the 
synchronous  averages  of  the  tail  rotor  gearbox  output  shaft  computed  using  these 
signals,  has  revealed  the  following: 

1.  There  is  significant  evidence  of  dynamic  twisting  of  the  tail  drive  shaft  during 
hover  and  tail  rotor  turns.  This  is  most  probably  caused  by  tail  rotor  torque 
variations  induced  by  pilot  control  pedal  movements. 

2.  The  level  of  this  twisting  is  significant  enough  to  affect  the  synchronous  averaging 
of  tail  rotor  gearbox  vibration  when  using  a  phase  reference  signal  derived  from 
the  main  rotor  gearbox. 

3.  During  forward  flight,  the  level  of  twisting  is  minimised,  and  both  phase 
reference  signals  perform  similarly  well  for  synchronous  averaging  purposes. 
However,  in  the  S-70A-9,  it  was  discovered  that  undesirable  vibration  impulses 
were  induced  in  the  tail  rotor  gearbox  vibration  during  these  flight  conditions.  It 
was  speculated  that  these  impulses  were  caused  by  the  aerodynamic  interaction  of 
the  main  rotor  tip  vortices  on  the  tail  rotor  blades. 

4.  The  (relatively  simple)  mathematical  model  developed  to  examine  the  theoretical 
effect  of  phase  reference  jitter  on  synchronous  averaging  proved  to  be  difficult  to 
apply  to  the  data  because  of  the  lack  of  an  absolute  phase  reference,  and 
assumptions  made  in  the  model.  Some  qualitative  agreement  with  the  model  was 
observed  for  the  forward  flight  data,  but  no  agreement  was  observed  for  the 
hover  flight  data. 

From  these  results,  the  following  recommendations  are  made: 

1.  Wherever  possible,  a  tail  rotor  phase  reference  signal  should  be  used  in  preference 
to  a  main  rotor  phase  reference  signal  when  synchronously  averaging  tail  rotor 
gearbox  vibration. 

2.  In  the  S-70  series  of  helicopters,  tail  rotor  gearbox  vibration  should  only  be 
examined  during  hover,  or  very  low  forward  speeds,  to  eliminate  the  potential  for 
any  undesirable  vibration  impulses  to  mask  faults  within  the  gearbox. 
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Appendix  1: 

S-70A-9  Black  Hawk  Transmission 


Shaft 

RPM 

Gear 

Teeth 

Mesh  Freq  (Hz) 

Input  Module 

High  Speed  Input  Shaft 

20900 

Input  Pinion 

22 

7663 

Input  Gear  (Quill)  Shaft 

5748 

Input  Gear 

80 

7663 

Accessory  Drive  Take-off 

76 

7281 

Accessory  Drive  Shaft 

11806 

Accessory  Drive  Gear 

37 

7281 

Accessory  Module 

Generator  Drive  Shaft 

11806 

Hydraulic  Drive  Take-off 

56 

11019 

Hydraulic  Drive  Shaft 

7186 

Hydraulic  Drive  Gear 

92 

11019 

Main  Module 

Combining  Pinion  (Quill)  Shaft 

5748 

Combining  Pinion 

21 

2012 

Crownwheel  Shaft 

1207 

Lower  Crownwheel  Gear 

100 

2012 

Upper  Crownwheel  Gear 

75 

1509 

Sun  Gear 

62 

980 
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Appendix  2: 

MATLAB  M  Files 


ATTEN.M 

function  atten (bw, sigma, conf,n) 

%atten (bw, sigma, conf ,n) 

%Plots  the  amplitude  response  of  a  signal  average  with  tacho  jitter. 

%bw  =  upper  frequency  limit  {shaft  orders) 

%sigma  =  max  standard  deviation  of  jitter  (degrees) 

%conf  =  confidence  level  [0,1) 

%n  =  number  of  averages 
%DMB  29/10/97 
f  req= (0 :bw) ; 
s= { 0 : 0 . 1 : sigma) *pi/180; 
conlim=mychi2conf (conf ,n) ; 
elf ; 

set (gca, ' colororder 1 , [0  0  0]); 

set (gca, 1 linestyleorder ' , 1 - 1  —  1 I : 1 ) ; 

set (gca, 'box1 , 'on'); 

set (gca, 1 ylim' , [0  1.05]); 

for  i=l : length ( s ) , 

t=-0 . 5*s (i) A2*f req. A2 ; 
a (i, : ) =exp (t) ; 

b(i, : )=exp(conlim(l) A2*t) ;  %LCL 
c (i, : ) =exp (conlim(2) A2*t) ;  %UCL 

11  (i,  : ) -sprintf  ( *  Wsigma  =  %1 .  If  Wcirc ' ,  s  (i)  *180/pi)  ; 

hp  (i)  =patch  {  [f  req  f  liplr  (freq)  ]  ,  [b  (i,  : )  ,  fliplr  (c(i,  :))]#-  [0.85  0.85 

0.85]); 
end 

hold  on 

h=plot (freq, a) ; 
hold  off 

xlabel (' Frequency  (shaft  orders)'); 
ylabel ( 'Amplitude ' ) ; 

title (sprintf ( 'Amplitude  Response  (%i%%  Conf.,  N  -  %i) ', conf *100 , n) )  ; 
legend (h, 11,3); 

SPEED.M 

function  speed ( tail, main, tl, t2,  yl,y2)  ; 

%  Calculate  and  plot  tacho  speed  variations 
fs=4000000/95;  %total  sample  rate 
fs=fs/3;  %sample  rate  per  channel 

time  =  (0 : length (tail) -1) ' /fs;  %create  time  vector 
main  =  main  *  1107225/21976768;  %  convert  to  tail  rotor  angle 
p=polyfit (time, tail, 1) ;  %linear  regression  of  tail  tacho 
y=polyval (p, time) ;  %evaluate  linear  fit 

vt  ~  (tail  -  y)*360;  %  subtract  linear  fit  and  convert  to  degrees 
vm  =  (main  ~y)*360;  %  subtract  linear  fit  and  convert  to  degrees 
vt  =  decimate (vt -mean (vt) , 100) ;  %subtract  mean  and  decimate 
vm  =  decimate (vm-mean (vm) , 100) ;  %subtract  mean  and  decimate 
time  =  ( 0 : length (vt ) -1 ) f*100/fs; 
figure (gcf ) ; elf; 

set (gcf , ' default axes font size ' , 8 ) ; 

set (gcf, 1 def aulttextfontsize ' ,  8)  ; 

set (gcf, ' defaultaxestickdirmode ' , 'manual ' ) ; 

set (gcf,  ' defaultaxestickdir ' , 'out' )  ; 

set (gcf , ' default axesbox ' , 'on'); 

set (gcf, 'position' , [100, 350, 680, 350]); 

subplot (1,2,1); 

hl=plot (time, vt , 1 k- 1 , time, vm, ' k — ' ) ; 
set (gca, 'XLim' , [0  max (time)]) 
set (gca, 1 YLim' , [-60  60] ) 
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xlabel ( 1  Time  (s )  ’ ) ; 

ylabel ( ’Angle  (deg) 1 ) ; 

title ( 1  Lead/Lag  of  Tail  Rotor  Angle1); 

legend (hi, *  Photocell  * ,  *  Generator  ’ ,  0)  ; 

xl  -  round (length (vt) *tl/max (time) ) ; 

x2  =  round (length (vt) *t2/max (time) ) ; 

subplot (1,2,2) ; 

h2=plot (time (xl : x2 ) , vt (xl : x2) , 1 k- 1 , time (xl : x2) , vm (xl : x2 ) , 1 k~- * ) ; 
set (gca,  ’ XLim' ,  [time (xl)  time  (x2) ] ) 
set (gca, 1 YLim' , [yl  y2] ) 

title (' Lead/Lag  of  Tail  Rotor  Phase  Angle1); 

xlabel ( 1  Time  (s) 1 ) ; 

ylabel ( ’Angle  (deg) ’ ) ; 

legend (h2, ’Photocell’ , ’Generator’ , 0) ; 

DPHASE.M 

function  delta  =  dphase (main, tail) ; 
main  =  main  *  1107225/21976768; 
delta  =  (tail-main) *360; 
delta  =  delta  -  mean (delta); 

PLOTD.M 

function  plotd (data, r) ; 
d=decimate (data, r)  ; 

time  =  (0 : length (d) -1) . ' *r*95*3/4E6; 

plot (time,d, '  k-  ’ )  ; 

set (gca, ’ ylim’ , [-1 . 2, 1 . 2) ) ; 

set (gca, ’xlim’ , [0, time (length (time) ) ] ) ; 

xlabel ( ’Time  (s) ' )  ; 

ylabel ( ’Angle  (deg)’); 

title (’Phase  Difference’); 

PDF.M 

function  pdf (delta, xmin, xmax, ymax, nbins) ; 

%  pdf (delta, xmin, xmax, ymax, nbins) 

%  Compute  probability  density  function  from  xmin  to  xmax  using  nbins. 
mu  =  mean (delta) ; 
sigma  =  std (delta); 

%  specify  mid  points  of  bins  for  hist  function 
dx  =  (xmax-xmin) /nbins; 
x  =  ( ( 0 : nbins-1 ) +0 . 5) *dx  +  xmin; 

%  compute  and  plot  histogram 
n  ~  hist (delta, x) ; 
p  =  n/sum (n) /dx; 
hs  =  stairs (x-0 . 5*dx, p, ’ k- ’)  ; 
set (hs, ’ linewidth’ ,  1)  ; 
haxes  =  gca; 

set (haxes, ’XLim’ , [xmin  xmax] ) 
set (haxes, 'YLim* , [0  ymax]) 

%  Label  graph 

set (gcf, 'position* , [10  320  500  400]) 
title ('PDF  of  Phase  Difference’); 
xlabel ( ’Angle  (deg)’); 
ylabel (’ Probability  Density’); 

%  Print  standard  deviation 
xtxt  =  xmin  +  0 . 95* (xmax-xmin) ; 
ytxt  =  0.9  *  ymax; 

str  =  sprintf  ( '  Wsigma  =  %6 . 3f\\circ ’,  sigma)  ; 
hsigma  =  text (xtxt , ytxt , str) ; 
set (hsigma, ' VerticalAlignment ’ , ’Top ' ) ; 
set (hsigma, ’HorizontalAlignment ’ , ’ right ’ )  ; 

%  Print  max  deviation 

str  =  sprintf  (’max  =  %  6 . 3f  Wcirc  ’  ,max  (delta)  )  ; 
ytxt  =0.8*  ymax; 


42 


DSTO-TR-0739 


hmu  =  text (xtxt, ytxt, str) ; 

set (hmu, 1 VerticalAlignment ' , 'Bottom1 )  ; 

set (hmu, 'HorizontalAlignment ' ,  ' right ' )  ; 

%  Print  min  deviation 

str  =  sprintf('min  =  %6. 3f\\circ min (delta) ) ; 
hsigma  =  text (xtxt , ytxt , str) ; 
set (hsigma, 'VerticalAlignment ' , 'Top' ) ; 
set (hsigma, 'HorizontalAlignment' , 'right' ) ; 

%  compute  and  plot  normal  distribution 
xx  =  xmin  +  (0 : 266) /266* (xmax-xmin) ; 
pn  -  1/sigma/sqrt (2*pi) *exp (-xx. A2/2/sigma. A2) ; 
hold  on; 

hpn  =  plot (xx,pn, ' k: * ) ; 
set (hpn, ' linewidth' ,0.5); 
legend (hpn, 'Normal  Distribution' , 2) 
hold  off; 

PSHIFT.M 

function  y  =  pshift  (x,  p) ; 

%  y  *  pshift  (x,  p) 

%  y  =  x  shifted  to  the  right  by  p  points  via  fft(x) 

%  x  =  signal 

%  p  =  number  of  points  to  shift 
fx  =  f ftshift (f ft (x) ) ; 

fz  =  exp(-j*2*pi* (-length(x) /2:length(x) /2-1) *p/length(x) ) ; 
fy  =  fftshift(fx  . *  fz)  ; 
y  =  if ft (fy) ; 

XCORR.M 

function  z  =  xcorr(sl,s2) 

%  z  ®  xcorr(sl,s2) 

%  z  =  circular  cross  correlation  of  si  and  s2  computed  via  fft 
%  si  =  signal  1 

%  s2  =  signal  2  (must  be  the  same  length  as  si) 

fl  -  fft (si) ; 

f2  -  fft  (s2) ; 

xspec  =  fl  . *  conj  (f2)  ; 

z  -  if ft (xspec) ; 

PSYNC1.M 

function  xcmax  =  psyncl (phi , si , s2 ) 

%  xcmax  =  psyncl (phi) 

%  xcmax  =  max  cross  correlation  of  signals  si  and  s2 
%  phi  =  number  of  points  to  shift  s2  to  the  right 
%  si  -  signal  1 
%  s2  =  signal  2 

%  index  =  point  of  max  correlation 
global  index 

[xcmax, index]  =  max (real (xcorr (si, real (pshift (s2, phi) )))) ; 
xcmax  =  -xcmax; 

COMPSIG2.M 

function  [phi, ds ] -comps ig2 (si, s2) ; 

%  [phi, ds] =compsig (si, s2) 

%  Compares  si  to  s2  by  finding  the  phase  difference  and  removing  it  from  s2. 
%  Plots  signal  and  power  spectra  of  si,  s2,  and  the  difference  signal. 

%  phi  =  number  of  points  that  s2  lags  si 
%  ds  =  difference  signal 
%  si  =  signal  1 

%  s2  =  signal  2  (must  be  same  length  as  si) 

%  index  =  point  of  max  correlation  computed  by  'psyncl' 

%  DMB  27/10/97 
global  index 

%  use  fmin  to  find  the  phase  shift  that  produces  the  max  correlation 
phi  =  fmin ( 'psyncl ' , 0, length (s2) -1, f opt ions, si, s2 ) ; 
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phi  =  mod (phi  +  index  -  1, length (s2) ) ; 

%  compute  difference  signal 
ds  -  si  -  real (pshift (s2, phi) ) ; 

%  plot  si 
figure (gcf ) ; elf ; 

set (gcf, ’position1 , [50, 60, 480, 600] ) ; 

set (gcf, 'defaultaxesfontsize* , 8) ; 

set (gcf, ' defaulttextfontsize  ' ,  8)  ; 

set (gcf, 1 defaultaxestickdirmode 1 , ’manual ’ ) ; 

set (gcf, ’ defaultaxestickdir ' , ’out’ ) ; 

set (gcf, ' default axesbox ’ , ’on’ ) ; 

x  -  360* (0 : length (si) -1) /length (si) ; 

f  -  (0 : length (si) /2-1) ; 

subplot (3,2,1) ,plot (x, si,  1  k* )  ; 

set (gca, ’xlim’ , [0,360]), set (gca, ’ylirn’ , [-15,15]); 

title ( ’ Synchronous  Average ’ ) , ylabel ( * g’ ) , text (-100,15, ’  (a) '); 

str  =  sprintf  ( *  Wsigma  =  %6.3f  gv , std(sl) ) ; 

ht  =  text (360, 15, str) ; 

set (ht, ’ VerticalAlignment  * , ’Top  * ) ; 

set (ht, 'HorizontalAlignment ’ , ' Right ’ )  ; 

fsl  =  f ft (si) ; 

subplot  (3,2,2); 

plot (f (1 : 400) , abs (fsl(l:400)) *sqrt (2) /length (fsl) , ’ k’ ) ; 
set (gca, ’xlim' , [0, 400] ), set (gca, ’ylirn* , [0,6]); 

%plot (f , abs (fs (1 : length (fs) /2) ) *sqrt (2) /length (fs) , ’ k’ ) ; 

%set (gca, ’xlim' , [0, length (fs) /2-1] ) , set (gca, ’ylirn' , [0,6]); 
title ( ’ Spectrum’ ) , ylabel ( ’ g  (rms) ’ ) ; 
for  1=1:5, 

str  =  sprintf (' %6 . 3f’ , abs (fsl (i*53+l )) *sqrt (2) /length (fsl) ) ; 
ht  =  text (f (i*53+l) , abs (fsl (i*53+l) ) *sqrt (2) /length (fsl) , str) 
set (ht , 'VerticalAlignment ' , ’Middle ’ ) ; 
set (ht, 'HorizontalAlignment ' , ' Left ’ )  ; 
set (ht , ’ Rotation ' , 90) ; 

end 

%  plot  s2 

subplot (3,2,3), plot (x, s2, ’ k’ ) ; 

set (gca, ' xlim' , [0,360]), set (gca, ’ylirn’ , [-15,15]) ; 

ylabel  Cg’)  ,  text  (-100, 15,  '(b)')  ; 

str  =  sprintf  (' Wsigma  =  %6.3f  g',std(s2)); 

ht  =  text (360, 15, str) ; 

set (ht , 'VerticalAlignment ' , ’Top’ ) ; 

set (ht, ' HorizontalAlignment ’ , ’ Right  * ) ; 

str  =  sprintf (’ Wphi  =  %6. 3f\\circ’ ,phi*360/length (s2) ) ; 

ht  *  text (360, 15, str) ; 

set (ht , 'VerticalAlignment ' , ’ Bottom’ ) ; 

set (ht , ’ HorizontalAlignment ’ , ’ Right ' ) ; 

f s2  =  f ft (s2 ) ; 

subplot (3,2,4) ; 

plot  (f (1:400) , abs (fs2 (1 : 4  00) ) *sqrt (2) /length (fs2)  ,  ’kf )  ; 
set (gca, 'xlim' , [0,400]) , set (gca,  ’ylirn’ ,  [0,6]); 

%plot (f , abs (fs (1 : length (fs) /2) ) *sqrt (2) /length (fs) , ’ k' ) ; 

%set (gca, 'xlim* , [0, length (fs) /2-1] ) , set (gca, ’ylirn’ , [0,6]); 
ylabel ( 'g  (rms) ' ) ; 
for  i=l:5, 

str  =  sprintf (' %6 . 3f' , abs (fs2 (i*53+l) ) *sqrt (2) /length (fs2)  )  ; 
ht  =  text (f (i*53+l) , abs (fs2 (i*53+l) ) *sqrt (2) /length (fs2)  ,  str) 
set (ht, ’VerticalAlignment' , ’Middle’ ) ; 
set (ht ,  'HorizontalAlignment ’ , *  Left 1 ) ; 
set (ht, ’Rotation’ , 90) ; 

end 

%  plot  difference  signal 
subplot (3,2,5) , plot (x, ds,  ' k’ )  ; 

set (gca, ’xlim’ , [0,360]), set (gca, 'ylirn* , [-5,5]); 

xlabel ( *  Phase  (degrees) ’ ) , ylabel ( ’g* ) , text (-100,5, ' (c)  ’); 

str  =  sprintf  (’ Wsigma  =  %6.3f  g',std(ds)); 
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ht  “  text (360, 5, str) ; 

set (ht, fVerticalAlignment 1 , 'Top1 ) ; 

set (ht, 'Horizontal Alignment' , 'Right' )  ; 

fs  =  f ft  (ds) ; 

subplot (3,2,6); 

%plot (f (1:400) ,abs (f si (1:400) ) . /abs (fs2 (1 : 400) )  ,  ' k' )  ; 
set (gca, 'xlim* , [0,400]), set (gca, 'ylim' ,[0,2.5]); 

%plot (f , abs (fs (1 : length (f s) / 2) ) *sqrt (2) /length(fs) , 'k') ; 

%set (gca, 'xlim' , [0, length (fs) /2-1] ) ,set (gca, 'ylim* ,[0,1.2]); 
xlabel (' Frequency  (orders) 1 ), ylabel ( 'Amplitude  Ratio'); 
for  i*=l:5, 

hl=line ( [f (i*53+l) ,f (i*53+l) ] , [0,abs (fsl (i*53+l) ) /abs (fs2 (i*53+l) )  ] )  ; 
set (hi, 1  color' ,[0,0,0]); 

str  =  sprint f ( ' %6 . 2f ' , abs (fsl (i*53+l) ) /abs (fs2 (i*53+l) ) ) ; 
ht  -  text (f (i*53+l) , abs (fsl  (i*53+l) ) /abs (fs2 (i*53+l) ) , str) ; 
set (ht, 'Vertical Alignment ' , 'Middle ' ) ; 
set (ht, ' HorizontalAlignment ' , 'Left' ) ; 
set (ht, 'Rotation' , 90) ; 

end 
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19.  ABSTRACT 

Synchronous  averaging  requires  an  accurate  phase  reference  (tachometer)  signal.  In  helicopter 
transmissions,  such  a  signal  can  usually  be  obtained  or  derived  from  an  engine  or  main  rotor  gearbox 
accessory,  but  is  not  as  readily  available  for  the  tail  rotor  gearbox.  This  report  examines  whether  a 
separate  tail  rotor  gearbox  phase  reference  signal  is  necessary  by  investigating  the  relative  jitter,  due  to 
dynamic  tail  drive  shaft  twist,  between  phase  reference  signals  obtained  from  the  main  and  tail  rotor 
gearboxes  in  a  S-70A-9  Black  Hawk  helicopter. 
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